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We present a quantized hydrodynamic theory and its applications of one-dimensional hard- 
core bosons in a harmonic trap. Quantizing the Hamiltonian of a trapped hard-core bosons and 
diagonalize it in terms of the phase and density fluctuations associated with the Bose field operator. 
As an applications, we calculate discrete energy spectrums, dynamic structure factor, momentum 
transferred by a two-photon Bragg pulse, single-particle density matrix, momentum distribution, 
and two-particle correlation function. The dynamic structure factor has multiple peaks due to the 
discrete nature of the eigen-modes which can be observed by a two-photon Bragg pulse with long 
duration. The coherence length at zero temperature is very small due to the "fermionization" of 
the system. We also compute the momentum distributions which has oscillatory behavior at large 
momentum. The pair distribution function shows that there are many deep valleys at various relative 
separations which implies shell structure due to the Pauli blocking in real space. 
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I. INTRODUCTION 



tion V, we present a brief summary of our work. 
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Recently, there have been renewed interest in the study 
of many-body effects in a one-dimensional(lD) quantum 
fluid. There is a realization of quasi-lD degenerate Bose 
gas in a harmonic trap [1] by freezing the radial mo- 
tion of atoms in a cylindrical trap to zero point oscilla- 
tions. Depending on the strength of the two-body inter- 
action, the trapped bosons are in two regimes: mean-field 
regime when the particles are weakly interacting which 
can be described by the Gross-Pitaevskii equation, and 
hard-core bosons (commonly known as Tonks-Girardeau 
gas) when the interaction is strong enough which can 
be described by the Lieb-Liniger model of ^-interacting 
ID bosons [2]. Most of the study concentrate on the 
phase fluctuating degenerate bosons in the mean-field 
regime [3-8] . With the possibility of creation of the hard- 
core bosonic system in a harmonic trap, there is a grow- 
ing interest to study theoretically this strongly interact- 
ing regime [9-18]. In this regime, the system acquires 
fcrmionic properties in a sense that there exists an exact 
mapping between the density of a strongly interacting 
bosons and the density of an ideal Fermi gas. The en- 
ergy per particle of hard-core bosons coincides with the 
Fermi energy of an ideal Fermi gas: E/N — (irhn) 2 / (6m) 
[19]. This is a manifestation of the Fermi-Bose duality in 
ID systems. 

This paper is organised as follows. In section II, we 
write down an effective hydrodynamic Hamiltonian of 
a trapped hard-core bosons and diagonalize it in terms 
of the phase and density fluctuations. In section III, 
we calculate the dynamic structure factor and the time 
evolution of the momentum transferred due to the two- 
photon Bragg pulse. In section IV, we calculate the 
single-particle correlation function, momentum distribu- 
tions and the two-particle correlation functions. In sec- 



II. QUANTIZATION OF ID HARD-CORE 
BOSONS 

We consider N interacting bosons in a cigar-shaped 
trap (u) z <C luq) a t very low temperature (T ~ 0). 
If temperature and the interaction energy per particle 
does not exceed the transverse level spacing hu>o, then 
the elongated cigar-shaped trapped system becomes ef- 
fectively one-dimensional. The two-body interaction in 
the axial direction is approximated by an effective two- 
body potential [17], U 1D = giS(z), where g\ = 
is an effective one-dimensional coupling constant and 

a 2 

dm = — -^-[1 — 1-4603 JS ] is the one-dimensional scat- 
tering length, where a is the scattering length in 3D sys- 
tem and an — w— ^— is the transverse oscillator length. 

The effective one-dimensional Hamiltonian for trapped 
atoms is 



where 



Ha 



H = H 



' 2m dzf 
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i=l 



N N 



S(zi - z 3 ) 



(1) 



(2) 



is known as for a ID ^-interacting Bose gas. This Hamil- 
tonian can be diagonalized via Bethe ansatz and the 
equation of state is known exactly for all densities and 
temperatures. The gas parameter is defined as the prod- 
uct of the density (n) with the effective scattering length 
aio, i.e n|ai£)|. At zero temperature, the interaction 



energy per particle is e(n) 



1 7T Tl 



i 2 when 



0. 



The low-density limit corresponds to the case of infinitely 
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strong interactions. We are interested to study the Bose 
systems with infinitely strong interactions at zero tem- 
perature by using the quantized hydrodynamic approxi- 
mation. Therefore, in the dilute regime, the Hamiltonian 
for the hard-core bosons confined in a harmonic trap at 
zero temperature can be written as, 

H' = J dz**{z,t){-^^ + ±rnu>y-ri*(z,t) 



9i 



dz|$(z,i)| 6 , 



(3) 



where g\ = (jT 2 h 2 )/(2m). The Thomas- Fermi density 
profile of the trapped hard-core bosons is different from 
bosons in the mean-field regime and it is 



n (z) = 




(4) 



where Zo = y/2Na z is the Thomas- Fermi length of the 
system. The chemical potential is fi = NTilo z . Also, one 
can define the Fermi wave- vector is kf — y/2N /a z — l/l c , 
where l c is the correlation length of this system. We ex- 
tend the quantized hydrodynamic theory developed by 
Wu and Griffin [6,20] at T = for the bosons in the 
mean-field regime to the hard-core bosons in a harmonic 
trap. It should be mentioned that the hydrodynamic 
theory of strongly interacting system is valid only when 
\z\ > l c or kl c < 1. The Bose order parameter can be 
written as ®(z,t) = \/n(z, t)e^ z '^ . The fluctuations 
in the density and phase about their equilibrium are 
n(z,t) — no(z) + Sn(z,t) and <j>{z,t) — (j>o(z) + S(j){z,t). 
By keeping up to second-order in the fluctuations, we 
obtain the effective hydrodynamic Hamiltonian at zero 
temperature, 



H = Hq 
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dz [mno(z)Sv 2 (z, t) + gin (z)Sn 2 (z, t)] 



where Sv z (z, t) = ^j^S(f>(z,t) is the superfluid velocity. 
This Hamiltonian is similar to the ID Hamiltonian de- 
rived by Haldane [21]. This Hamiltonian is quadratic in 
the density and the phase fluctuation operators, then we 
can diagonalizc it using the canonical transformations: 

5n(z, i) = £ [a j f j (z)e- i ^ t a j + H.c] , (5) 



and 



54>(z,t) = £ [b^zy-W&i + H.c] , (6) 



where fj{z) 



, The operators d, and a\ de- 



stroy and create excitations with energy Tiujj and satisfy 
the commutation relations [aj,&L] = [6(j,&j] = 



and [&],&]] = 0. Also, [Sn{z,t),S$(z',t)] = iS(z - z'). 

The constant terms are a,j — i\J 4 g ^no( W) anc ^ = 

\J 9l h^ ■ e q ua ti° ns f° r the density and phase fluc- 
tuations can be obtained by using the Heisenberg equa- 
tion of motion: 



d 2 6h(z,t) 

dt 2 



d_ 

dz 



—n (z)-^-[n (z)Sh(z,t)] 
m dz 



and 



d 2 5cj){z,t) 
dt 2 



—n (z)-^-[n (z)S4>(z,t)} 
m dz 



The equation for the eigenfunctions ipj{z) is 
[(1_52) £"^ + ^ (Z) = ' 



(7) 



(8) 



(9) 



where z = z/Zq. Therefore, Eq. (9) becomes Chebyshev 
equation with the eigenfrequencies given by ujj — joj z and 
the corresponding normalized eigenfunctions are ipj(z) — 

\J~-^Tj (z I Z ) , where Tj(z/Z ) is the first-order Cheby- 
shev polynomial in z and So — 1, otherwise Sj = 2. The 
monopole mode corresponds to the j = 2 mode which is 
independent of the interaction strength [22] due to the 
conformal symmetry [23] of the Hamiltonian (see Eq. 
(3)). Note that the eigenfunctions of the density and 
phase fluctuations are different in the hard-core bosons 
regime, in contrast to the case of bosons in the mean-field 
regime. 



III. STRUCTURE FACTORS AND BRAGG 
SPECTROSCOPY 

Dynamic Structure Factor. Consider a low-intensity 
off-resonant inelastic light scatters with momentum 
transfer fik and energy transfer hw to the target (hard- 
core bosons). If the external light couples weakly to the 
number density of the target, the differential cross section 
is proportional to the dynamic structure factor S(k,u>), 
which is obtained from the Fourier transform of the time- 
dependent density-density correlation functions, 



S(k,w) 



J dt J dze l{ut ' kz) < Sh(z,t)Sh(0,0) > . (10) 



It is the density fluctuation spectrum that can be mea- 
sured in the two-photon Bragg spectroscopy. The dy- 
namic structure factor of this system can be written as, 



^4^0(0)'^ 
x [NjS(u; + ujj) + (1 + Nj)S(lo - uj)], 



(11) 
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where Nj = [exp((3hujj) — is the thermal Bosc- 

Einstein function and Fj(k) = J dze~ lqz fj(z) is the spa- 
tial Fourier transformation of fj(z). It can be easily 
shown that 



\F 3 (~k)\ 2 = Z n 2 \J^N~k)\ 2 



where the dimensionless variable is k = k/kf. 
can write the dynamic structure factor at T = 
S(k,u>) — V ■ Sj(k)S(uj — u)j)], where 



S 3 {k)=j\J 3 {2Nkt 



(12) 



We 
as 



(13) 



These functions determine the weight of the light- 
scattering cross-section in S(k, u>) of the corresponding 
collective modes of energy Wj. In Fig.l we plot Sj(k) 
as a function of the dimensionless wave vector k = k/kf 
for the excitations j = 1, 2, 3, 4, and 5. Fig.l shows that 
for a given k, how many modes significantly contribute 
to S(k,uj). Note that j = mode do not contribute 
to S(k,uj). It is clear from the Fig.l that the strongest 
weights for these collective modes appear for k > 0.02, 
it implies that the momentum transfer k in a light scat- 
tering experiments should be k > 0.18OJ 1 for N = 40 in 
order to pick up the strong spectral weight from the low- 
energy collective modes. We also notice that the number 
of modes that contribute to the S(k,oj) increases with 
the chemical potential. 
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FIG. 2. Plot of the dynamic structure factor S(k, ui) vs. 
u/uj z at T = 0. 

As one can see from the Fig. 2, the dynamic structure 
factors has multiple peaks. This multiple peaks are due 
to the underlying discrete spectrum. Note that the en- 
ergy transfer lu is much smaller than the gap between 
the chemical potential \x and the transverse excited state 
(~ ojo) and hence the u) can not excite the radial modes 
and its contribution to the dynamic structure factor is 
zero. This behavior of the multiple peaks can be re- 
solved in two-photon Bragg spectroscopy, as shown by 
Steinhauer et al. [24]. 

Bragg Spectroscopy: The observable in the Bragg scat- 
tering experiments is the momentum transferred to the 
system. When the system is subjected to a time- 
dependent Bragg pulse, the additional interaction term 
appears in the Hamiltonian which is given by [25], 



Hi(t)= dzip^(z,t)[V B (t)cos(kz-u}t)]ip(z,t). (14) 



Here, we supposed that the Bragg pulse is switched on 
at time t = and k is along the z-direction. 

The momentum transfer from the optical potential is 
obtained in Ref. [8] is given by 



P z (t)=^k<a](t)a j (t)>= ( 



v B {t) 
2h 



FIG. 1. Plots of Sj(k) as a function of the dimensionless 
wave vector k for iV = 40. 



x F 3 [{uj - u), t] - Fj [{uij +u),t] 
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(15) 



For k — > 0, the leading term arises from the Kohn mode 
j = 1, with Si ~ fc 2 ; the next contributions arise from 
the terms with j — 2, 3. In Fig. 2 we plot the dynamic 
structure factor S(k,ui) for the momentum transfer cor- 
responds to k = O.lfc/. For finite-energy resolution we 
have replaced the delta function by the Lorentzian with 
a width of T = 0.12w z . 



wherei^±u;),i] = ( ""fetT/f ) 2 ^Sj(k) is given 
in Eq. (13). For large t, S(k,uj) - P z (t). We plot 
P z {t) for various time in Fig. 3 and shows that the mul- 
tipeak spectrum in S(k,u>) can be resolved only when 
the duration of the Bragg pulse is t » 2tt/uj z . When 
t < 2tt/oj z , P z {t) reflects the dynamic structure factor 
calculated from the local density approximation. 
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FIG. 3. (Color online). Plots of the momentum transferred 
Pz(i) vs. lu at T = when momentum transfer is k = O.lfc/ 
for various time ui z t = 4.0 (dashed), u z t — 10.0 (dotted) and 
uj z t = 19.0 (solid). Also, we have assumed Vb = 0.3hui z 



Static Structure Factor. For completeness, we also cal- 
culate the static structure factor which can be obtained 
from the Fourier transform of the equal-time density- 
density correlation function 

S(k)= [ dze lk( - z - z,) (5n(z)5h(z')) 



(16) 



with the long-wavelength limit S{k) ~ k 2 coth(!-j±). 

IV. CORRELATION FUNCTIONS AND 
MOMENTUM DISTRIBUTION 

The equal-time single- particle density matrix is defined 



D 1 (z,z') = (ft(z,t)$(z',t)). 



(17) 



The single-particle density matrix by taking care of the 
density and the phase fluctuations upto quartic term is 
given by 



where 



F 2 [z,z']= J Y,-[T 1 {~z)-T ] {~z')] 2 coth 



3=0 



, (18) 



and Fi[z, z'\ = Fftz, z'] - 3(F 2 [z, z'}) 2 , where 



Fi[z, z') = J2- ~ T ^')] 4 (! + 2N J + 2N D- 

j=o 3 



and 

F d [z,z'] 



^ AN 2 

3=0 



Tj(z) T 3 (z') 



1-z 



?2 



1-Z 



-1 2 

coth 


f3hujj 
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We calculate all the summation within the phonon 
regime, hujj < fi and the upper cut-off limit, j max , is 
obtained from the relation, [i = 7ia; Jmax i.e. j max = N. 
The normalized one-body density matrix or the phase 
correlation function, Ci(z), is denned as 



(ftWfc')) -j[F 2 -% + %] (1Q) 

y/n (z)n (z') 



It is known that the strong two-body interaction induces 
large phase fluctuations, and therefore we must consider 
the higher order term in the phase fluctuations to pro- 
vide more accurate coherence function. In fact, we have 
checked that the quartic phase fluctuation term (F4) re- 
duces the phase coherence function significantly (15-20 
percentage). We plot the phase coherence C\{z) vs. the 
separation z in Fig. 4. 
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FIG. 4. (Color online). Plots of the normalized one-body 
density matrix Ci(z) vs. the the separation 5 for various 
number of particles at T = 0. 

At T = 0, the coherence function goes to zero even 
at small separations due to the "fermionization" of the 
strongly interacting bosons. It has oscillatory behavior 
at large z which is also reflected in the momentum dis- 
tributions (see Fig. 5). Note that the coherence function 
decreases as the number of particles increases [15]. 

The momentum distribution can be obtained from 
Fourier transformation of the one-particle density ma- 
trix: 



n{k) 



° dz f ° dz'Diiz^'y^-*"). (20) 

— Zn J— Zn 



Fig. 5 shows that the momentum distribution has a os- 
cillatory behavior at large momenta which can be ob- 
served in the Bragg spectroscopy [26]. Note that this 
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hydrodynamic approximation is valid only when kl c < 1 
i.e. kZo < 2N. The momentum distributions for dif- 
ferent number of particles which is shown in Fig. 5 falls 
within the valid hydrodynamic regime. 
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FIG. 5. Plots of the zero temperature momentum distribu- 
tion n(fc)/n(0) vs. momentum k = kZo for N — 10 (solid) 
and N = 5 (dot-dashed). 

The width of the momentum distributions decreases 
with the particle numbers [13,15]. The oscillatory be- 
havior in the momentum distribution is also obtained by 
exact calculation in Ref. [15]. Fig. 4 shows that the single- 
particle correlation functions follows the power-law decay 
at large distances:Ci(z) <~ l/(y/z). Also, Fig. 5 shows the 
momentum distributions follows the power-law decay at 
large momentum: n(k) ~ l/(\/fc). These power-law de- 
cay in this system has been recently observed by using 
the exact quantum Monte Carlo techniques [16]. 

Similar to the single particle correlation function, we 
also calculate two-point correlation function at T = by 
taking care of the phase fluctuations upto quartic term 
which is given by, 

C 2 {z 1 ,z 2 ,z 3 ,z 4 ) = e ^ + s (21) 

where 

N 1 

G 2 [z t ] = -PM^i) + T J&) - T Ah) - Tj(zi)] 2 (22) 
i=i 3 

and 

N 1 

G ^ = E pFifr) + T J&) - T Ah) - T,(z 4 )] 4 

3=1 3 

- (G 2 N) 2 - (23) 

We plot the C 2 vs. the relative separation s for differ- 
ent values of mean-distances d in Fig. 6. Here, s and 
d (in units of Zq) are defined as z\ = (d + s)/2,z 2 = 
(-d-s)/2,z 3 = {-d + s)/2, and z A = (d-s)/2. 



C 2 (s) 
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FIG. 6. Plots of the zero temperature two-point correla- 
tion function C2(§) vs. the relative separations s = s/Zo for 
various mean-distances d 



The equal-times pair distribution functions shown in 
Fig. 6 contains information on the relative spatial distri- 
bution of pairs of hard-core bosons. We should men- 
tion that Fig. 6 is valid only when s > 1/(2N) due the 
hydrodynamic approximation. Actually, C 2 should van- 
ish with the relative distance s due to the hard-core na- 
ture of the particles. This hydrodynamic approximation 
is failed to describe this short-range correlations due to 
the strong interactions when s < 1/(2N). The effect 
of Pauli exclusion principle in real space is also observ- 
able as a deep valley at various relative distances s. It 
implies that the strongly interacting bosons in ID trap 
behaves as a fermions. This effect is absent in a quasi- 
1D Bose gas in the mean- field regime [27]. At large dis- 
tances the pair distribution function goes rapidly to zero 
and it decouples into the product of two-particle den- 
sity profiles. Recently, the two-point correlation function 
has been measured experimentally in an elongated de- 
generate Bose gases [28] and this method can be used to 
observe the valleys in the two-point correlation function 
of a trapped hard-core bosons. 



V. SUMMARY 



In this work, we have carried out a detail study of 
hard-core bosons in a harmonic trap by using the hy- 
drodynamic approximation at zero temperature. Par- 
ticularly, wc have calculated dynamic structure factor, 
single-particle correlation function, two-particle correla- 
tion functions and momentum distributions. We found 
many interesting behavior in the observables which are 
absent in the trapped bosons in the mean-field regime. 
All the results are obtained analytically. These results 
can be used to identify the hard-core bosons in a trap 



5 



from the bosons in the mean-field regime, and also im- 
portant with the progress of the experimental technique 
to achieve this strongly-interacting regime. 

I would like to thank M. A. Cazalilla for carefully read- 
ing the manuscript and valuable comments. 
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